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ABSTRACT 

A new algorithm for five-hole probe calibration and data re- 
duction using a non-nulling method is developed. The significant 
features of the algorithm are: 1) two components of the unit vector 
in the flow direction replace pitch and yaw angles as flow direc- 
tion variables, and 2) symmetry rules are developed that greatly 
simplify Taylor’s series representations of the calibration data. 
In data reduction, four pressure coefficients allow total pressure, 
static pressure, and flow direction to be calculated directly. The 
new algorithm’s simplicity permits an analytical treatment of the 
propagation of uncertainty in five-hole probe measurement The 
objectives of the uncertainty analysis are to quantify uncertainty 
of five-hole probe results (e.g., total pressure, static pressure, and 
flow direction) and determine the dependence of the result un- 
certainty on the uncertainty of all underlying experimental and 
calibration measurands. This study outlines a general procedure 
that other researchers may use to determine five-hole probe result 
uncertainty and provides guidance to improve measurement tech- 
nique. The new algorithm is applied to calibrate and reduce data 
from a rake of five -hole probes. Here, ten individual probes are 
mounted on a single probe shaft and used simultaneously. Use 
of this probe is made practical by the simplicity afforded by this 
algorithm. 
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Figure 1 A typical five-hole probe 


(j) = yaw angle 

ip = probe outer tube opening angle 

£ = rake probe circumferential position 

INTRODUCTION 

Figure 1 shows a typical five-hole probe and the coordinate 
system used in this analysis. The openings of the four outer tubes 
are inclined from the normal to their axes by angle while the 
center tube opening is normal to its axis. The pressure sensed 
at the opening of each tube is denoted p, . The subscript i refers 
to the probe opening, as defined in Figure 1. The mean velocity 
vector is V = f/i + Vj-f Wk and the unit vector in the mean flow 
directions is V /|V| = ui + vj + ink. Table 1 prescribes the unit 
vectors in terms of pitch and yaw angles in both the pitch-yaw 
and yaw -pitch systems of angle definitions. The sign convention 
was chosen so that (small) positive values of pitch and yaw angle 
provide positive v and w velocity components. When flow angles 
are required in this analysis the pitch-yaw system is used however 
the conclusions are valid for the yaw-pitch system. 

The ability to measure total and static pressure (hence veloc- 
ity magnitude) and flow direction with a five-hole probe is well 
established. Two different measurement techniques have been em- 
ployed. For the yaw-nulling method, the flow yaw angle is de- 
termined by the amount of rotation required until the pressures 
measured by tubes 2 and 4 are equal. Flow pitch angle is de- 
termined by the pressures of tubes 1, 3, and 5 and the results of 
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sin 6 cos <j> 


Table 1 Unit vectors 


pitch angle calibration. For the non-nulling method, the probe 
remains stationary and calibration is used to determine flow direc- 
tion. 1 In both methods, calibration is used to determine total and 
static pressure. Each method has advantages. Calibration and data 
reduction require less effort with the yaw-nulling method. When 
acquiring data in an experiment, the yaw-nulling method requires 
more time and more sophisticated actuation and data acquisition 
hardware than the non-nulling method. When the amount of data 
gathered in an experiment is large, the extra effort expended to 
calibrate the probe for the non-nulling method is more than com- 
pensated by the effort saved when acquiring data. This analysis 
is about the non-nulling method. 

Four pressure coefficients. Equations 1-4, are used in the cali- 
bration and data reduction procedure. Over a range of Reynolds 
numbers and subsonic Mach numbers the four pressure coefficients 
depend on flow direction only. Calibration is used to determine 
the relationships that exist between the four pressure coefficients 
and flow direction. In data reduction, Cp v and Cp w allow flow 
direction to be determined, while Cps and Cp a v g allow total and 
static pressure to be determined. 

Cp avg = (1) 

Po-P 

Cps = ELZl ( 2 ) 

P o - p 

C Pv = JLZIL. ( 3 ) 

Pb Pavg 

Cpw = (4) 

Pb Pavg 

There continues to be a great deal of interest regarding the need 
to quantify and report uncertainty in numerical and experimental 
fluid flow studies (Rood and Telionis, 1991 and Sutton, 1994). 
Most five-hole probe calibration and data reduction algorithms 
are not generally suited for an analytical treatment of uncertainty 
analysis. Uncertainty may be assessed by propagating measurand 
uncertainty through the data reduction programming by a numer- 
ical technique known as "jitter” analysis (Moffat, 1982). While 
the jitter technique allows the uncertainty of individual probes to 
be quantified, it provides little insight into basic processes un- 
derlying the uncertainty that analytical treatments yield. In the 
following section we develop an algorithm with a simplicity and 
accuracy that permits an analytical treatment of the propagation 
of uncertainty, revealing details about the relationships between 
probe geometry, calibration, data reduction and result uncertainty 
and providing guidance to improve measurement technique. 


1 The choice of pitch and yaw angle definitions depends on 
several factors. For the non-nulling technique the angle definition 
depends on the probe actuation hardware used in calibration. For 
the yaw-nulling technique the yaw-pitch angle definition must be 
used. 
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THE NEW ALGORITHM 

The discussion up to this point is classical. Excellent descrip- 
tions are given by Bryer and Pankhurst (1971) and Treaster and 
Yocum (1979). Although other authors use different names and 
definitions, common to all these methods is the use of four pres- 
sure coefficients, two coefficients to determine flow direction and 
two coefficients to determine total and static pressure. Our algo- 
rithm differs from previous research in our choice of flow direction 
variables and how the relationships between the four coefficients 
and flow direction are formulated. 

Traditionally, pitch and yaw angles are the variables used to 
express flow direction. Instead, we use components of the unit 
vector in the flow direction. Because the three components satisfy 
the criteria u 2 -1- v 2 + w 2 = 1, only two are required to uniquely 
specify flow direction (with the restriction that u > 0). We use 
v and w. Using two unit vector components instead of pitch and 
yaw angles has several advantages: 1) the calibration results do not 
depend on the choice of pitch-yaw or yaw-pitch angle definitions, 
2) probe tip symmetry is expressed easier, 3) the relationships 
between pressure coefficients and flow direction are simpler and, 
4) a step is removed from the data reduction procedure, where the 
desired answer is usually the velocity vector and not pitch and 
yaw angles. During calibration pitch and yaw angles are recorded 
and are later converted to unit vectors components according to 
the definitions in Table 1. 

Two general classes of schemes that have been reported to ex- 
press the relationships between the four pressure coefficients and 
flow direction are graphical methods (Dudzinski and Krause, 1969 
and Bryer and Pankhurst, 1971) and interpolative methods involv- 
ing spline fits of all the calibration data (Treaster and Yocum, 1979 
and Lagrani et ai. y 1989). Graphical methods are not suitable for 
large experimental data sets and the interpolative methods become 
cumbersome and slow when the amount of calibration data is large. 

We use Taylor’s series expansions in two variables, truncated 
to fifth order, to approximate the relationships between the four 
pressure coefficients and flow direction. What is unique in our 
approach is to identify and utilize the probe symmetry to reduce 
the number of required coefficients for all Taylor’s series from 
sixty to thirteen. For data reduction, flow direction as well as 
total and static pressure are obtained from the probe data by direct 
calculation. No piecewise interpolative or iterative procedures are 
required. 

An examination of Figure 1 will show that the five-hole probe 
possesses three types of geometric symmetry. They are: 1) sym- 
metry with respect to reflection across the xz-plane, 2) symmetry 
with respect to reflection across the xy-plane, and 3) symmetry 
with respect to 90° rotation about the x-axis. The symmetry of 
the probe has important consequences for the relationships be- 
tween the pressure coefficients and the flow direction variables. 
One result of symmetry is that within a range of flow directions 
about v = 0 and w = 0 there exists a one-to-one relationship be- 
tween the flow direction variables v and w and the flow direction 
coefficients Cp v and Cp w . In fact, the ability to use a five-hole 
probe to measure flow direction requires the existence of this one- 
to-one relationship. This allows either v and w or Cp v and Cp w 
to be used as independent variables to express flow direction. 


Since Cp av g and Cpb depend on v and w , the one-to-one 
relationship described allows Cp a vg and Cp$ to alternatively be 
considered functions of Cp v and Cp w . In Equations 5 and 6 the 
coefficients Cp avg and Cp$ are represented as a Taylor’s series 
expansions in the variables Cp v and Cp w ? 

oo oo 

Cp<xvg{Cpv > G’pty) = ^ ^ ^ ^ &ijCp v Cpw (5) 

i=0 j = 0 
oo oo 

Cps(Cp Vy Cp w ) = 'Y^2,bijCp' v Cp J w (6) 

»=0 j=0 

Symmetry relationships for Cp a vg and Cps are listed below. 
Proofs of the results stated below are included in the Appendix: 

Cpavg (Cp v , Cpxv ) = ^?avj( C?p v , Cpw ) 

Cpavg(Cpvi Cp w ) = Cp aV g(Cp v , Cpti;) 

CPavg{Cp v , Cpw) = Cp av g(C/p w , C p v ) 

Cp$ > Gpu; ) = Cp$( GpvjCpu/) 

Cps(Cp v y Cpw^) = G’ps(C , p v , Cp w ) 

Cpb (Gpv > Cp w ) = CpbiCpw) G*p v ) 

The symmetry conditions listed for the coefficients Cp a v g and 
Cpb greatly simplify their Taylor’s series representation. Also, 
subsonic flow aerodynamics requires Cps(0, 0) equal one. In the 
Appendix, the following results for the Taylor’s series are derived: 

dij = 0, j — 1)3)5,..., 

CLij =0, I = 1,3,5,..., 

dij = dji 

bij — 0, j = 1,3,5,..., 

bij — 0 , { = 1,3,5,..., 

bij — bji 

The simplified Taylor’s series for Cp a v g and Cps. truncated 
at fifth order and given by Equations 7 and 8, result in seven 
undetermined series coefficients, rather than thirty that would 
generally be required. The Taylor’s series coefficients have been 
renumbered in Equations 7 and 8 for clarity. The “hat” symbol 
that appears above the pressure coefficient is used to distinguish 
the Taylor’s series approximation from the pressure coefficient 
based on actual measurement. 

Cpavg =<*1 + 

a 2 {Cpl+Cpl)+ (7) 

a 3 {Cpl + Cpi,) + a+Cp 2 v Cp 2 w 

Cpb =14- 

bi{Cpl+Cpl)+ (8) 

b 2 {Cp* v + Cpl) + b 3 CplCpl 


2 Using Cp v and Cp w as independent variables in Equations 
5 and 6 allows Cp aV g and Cp$ to be calculated directly from Cp v 
and Cp w . If v and w had been used as independent variables in 
Equation 5 and 6, then an additional step would be required in 
data reduction when determining Cp Qvg and Cp$. 
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For flow direction, the unit vector components v and w are 
represented in Equations 9 and 10 by Taylor’s series expansions 
in the variables Cp v and Cp w ? 


v{Cp v ? Opw) — ^ ^ ^ CjjC'Px) Cpj, (9) 

i=0 j=0 

oo oo 

v)(Cp v , Cpw') — EE dijCpiCpi (10) 

»=o j=0 

In the Appendix the following probe symmetry results for v 
and w as functions of Cp v and Cp w are derived. 

v(Cp v ,Cp w ) = - v(-Cp v ,Cp w ) 

v{Cpv j (-'Ptu) = v(Cp v , C'pu/) 

w(Cp Vl Cp w ) = w(—Cp v ,Cp w ) 

w(Cp vy Cp w ) = —w(Cp v1 —Cp w ) 

v{Cpv > Cput) — w(Cpxjj , Cp v ) 

In the same manner used for Cp a v g and Cps, the results listed 
above greatly simplify the Taylor’s series representation of v and 
w. In the Appendix, the following results for the Taylor’s series 
coefficients are derived. 

Cij = 0} i = 0, 2, 4, . . . 

Cij = 0 , j = 1 , 3 , 5 ,... 

dij — 0, i = 1, 3, 5, . . . 

dij = 0, j = 0,2,4,... 

Cij = dji 

Equations 11 and 12 give the Taylor’s series approximation of v 
and w truncated at fifth order. The series coefficients have been 
renumbered in Equations 11 and 12 for clarity. The same series 
coefficients appear in both Taylor’s series expansions. Rather 
than requiring thirty coefficients, the truncated Taylor’s series 
only requires six. The “hat” symbol again is used to distinguish 
Taylor’s series approximations from actual values. 

v(Cp Vy Cp w ) =ciCpt,+ 

C2Cp Z v -f CzCpvCplj^- (11) 

c\Cp b v + c$Cp z v Cp 2 w + ceCpvCpw 

w(Cp Vf Cp w ) =ciCp w + 

czCp 2 v Cp w + C 2 Cp z w + (12) 

c$Cp\Cp w + csCplCplv -f c\Cp b w 

Note that all the desired simplified Taylor’s series, given by 
Equations 7, 8, 11, and 12 include terms involving products of 
Cp v and Cp w - In order to evaluate these coefficients requires 
taking calibration data over a range of both pitch and yaw angle 

3 Note that v and w are dependent variables in Equations 9 
and 10. In data reduction, this allows v and it; to be calculated 
directly from measured values of Cp v and Cp w . If Cp v and Cp w 
were expressed as functions of v and w, then a numerical solution 
scheme would be required to calculate values of v and w from 
measured values of Cp v and Cp w . 


values. In other words, it is not sufficient to simply vary pitch 

while yaw equals zero and vice-versa. 

The coefficients themselves are found by a least-squares proce- 
dure. This results in three systems of simultaneous linear equa- 
tions, the solutions of which are the Taylor’s series coefficients. 
The three systems of equations are given in Equations 13, 14 and 
15. The summation in these equations is over all data acquired 
during calibration. These equations are solved directly by matrix 
inversion. 
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(15) 


There are several additional details included in the calibration 
and data reduction algorithm. Figure 1 shows a Cartesian coor- 
dinate system relative to the probe tip. It is useful to think of 
a second Cartesian coordinate system relative to the probe shaft. 
The probe shaft coordinate system has its z-axis parallel to the 
probe shaft and its x-axis is the datum used for measuring probe 
shaft rotation. In practice, all measurements of the probe orien- 
tation are relative to the probe shaft coordinate system. Ideally, 
the probe tip is mounted so that the two coordinate systems have 
identical orientations. In reality this is highly unlikely. Four 
parameters are needed to account for misalignment between the 
probe tip and probe shaft. They are pitch and yaw offset and two 
tip rotation parameters. 

Pitch and yaw offset relate the difference between the x-axes of 
both coordinate systems. The value of pitch offset is determined 
by the probe geometry and should not change unless the probe is 
damaged. The yaw offset value is established whenever the probe 
is attached to the calibration rig or experimental rig. Thus, the 
yaw offset may change whenever the probe is moved. 

The tip rotation parameters relate the difference between the y- 
and z-axes of both coordinate systems when the x-axis difference 
has been eliminated. This is shown in Figure 2. To be rigorously 
correct, if the two tip rotation parameters are such that v\ ^ w \ , 
then the probe tip is skewed and the symmetry criterion that 
were used to develop the calibration procedure are not satisfied. 
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Figure 2 Probe tip rotation angles 

However, small amounts of skewness, which are unavoidable, 
result in only a small deviation from the symmetry criterion. 
Measuring and allowing for minor skewness in this manner can 
be viewed as a perturbation of the symmetry condition. The two 
tip rotation angles are also determined by the probe geometry and 
are unlikely to change unless the probe is damaged. 

Values of the four parameters are determined from calibration 
data. This is necessary because it is impractical to accurately de- 
termine these quantities from observations of the probe geometry 
alone. First, pitch and yaw angles recorded during calibration are 
converted to unit vector components. The pitch offset is deter- 
mined by evaluating the data recorded for zero yaw angle (hence 
v = 0) to find the value of w when Cp w = 0. This value, tu 0 , is 
the pitch offset A least square error procedure is used to deter- 
mine its numerical value. Likewise, the yaw offset is determined 
by evaluating the data recorded for zero pitch angle (hence w = 0) 
to find the value of v when Cp v = 0. Equations 16 and 17 trans- 
form unit vector components for pitch offset (Equation 16) and 
yaw offset (Equation 17). The primed unit vector components 
indicate values uncorrected for offset. 
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Determining the tip rotation angles is slightly more complicated. 
If there were no tip rotation, then the vector VCpt/(0,0) would 
point in the direction of the y-axis and the vector VCpu,(0,0) 
would point in the direction of the z-axis. The tip rotation 


parameters are a measure of the angles between the two gradient 
vectors and their corresponding axes. They are calculated using 
Equation 18. Values of the partial derivatives are determined 
with another least squares procedure. Equation 19 is the final 
relationship applied to transform unit vector components from the 
probe shaft coordinate system to the probe tip coordinate system. 
The primed unit vector components indicate values uncorrected 
for tip rotation. 



In summary, the steps used in the probe calibration algorithm 

are: 

1. Express the measured pressures as pressure coefficients 
Cpavg . Cps, Cp v and Cp w using Equations 1, 2, 3 and 
4. 

2. Convert the pitch and yaw angles recorded during calibration 
to unit vector components using the definitions from Table 1. 

3. Determine pitch and yaw offset w 0 and v 0 from the calibra- 
tion data and apply them to the unit vector component values 
using Equations 16 and 17. 

4. Determine the tip rotations parameters vi and w\ from the 
calibration data using Equation 18 and apply them to the unit 
vector component values using Equation 19. 

5. Calculate the Cp a vg(Cp Vy Cp w ) Taylor’s series coefficients 
using Equation 13. 

6. Calculate the Cps(Cp Vy Cp w ) Taylor’s series coefficients 
using Equation 14. 

7. Calculate the v(Cp Vl Cp w ) and w(Cp Vl Cp w ) Taylor’s se- 
ries coefficients using Equation 15. 

Once the probe is calibrated, the procedure for reducing ex- 
perimental data is straightforward. The steps that are applied are 

listed below: 

1. Determine the coefficients Cp v and Cp w from the probe 
pressure data using Equations 3 and 4. 

2. Substitute the values of Cp v and Cp w into Equations 11 
and 12 to calculate v and w. 

3. Substitute the values of Cp v and Cp w into Equations 7 and 
8 to calculate Cp aV g and Cp$. 

4. Calculate total and static pressure with Equation 20, which 
is simply a rearrangement of Equations 1 and 2 (with Cp avg 
and Cps substituted for Cp a vg and Cp$). 
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Figure 3 Pitch and yaw calibration results 


Figure 4 Cp a v 9 calibration results 


5. Convert the unit vector components from the probe tip co- 
ordinate system to the probe shaft coordinate system using, 
in order. Equations 21, 22 and 23. These equations are the 
inverse of Equations 19, 17 and 16. 
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( 21 ) 


( 22 ) 


(23) 


constant w nominally horizontal. An examination of Figure 3 will 
verify the validity of Equations 11 and 12. Average (RMS) error 
for both v and w is less than 0.0053 (or less than 0.3°). Likewise, 
contour plots of the Cp a vg and Cp$ coefficients as functions of 
independent variables Cp v and Cp w are shown in Figures 4 and 
5. Solid lines are the Taylor’s series approximations, broken lines 
are calibration data. Average error for both coefficients is less than 
1% of dynamic pressure. Again, the agreement is quite good. 

Probe Geometry and Response 

By making some plausible assumptions it is possible to re- 
late the leading order probe calibration coefficients for Cp a vg » 
v and w to probe geometry. The leading order term for Cp$ 
was determined earlier from aerodynamic considerations, that is 
Cps( 0,0) = 1. This analysis leads to the following results; 
a\ » cos 2 rl> and ci % ~ sin ^ • These results are approximate 
and are not intended to replace calibration, but they do provide 
some insight into the effect of probe geometry on probe response 
and uncertainty. 


It is worth reiterating the fact that neither the calibration or data 
reduction procedures requires iteration. All results are obtained 
by direct calculation. 

Experimental Verification of the Algorithm 

It remains to show that the simplified Taylor’s series, given 
by Equations 7, 8, 11 and 12, whose derivation was based on the 
symmetry of an ideal probe, can faithfully represent the calibration 
data of an actual five-hole probe. Flow direction calibration results 
for a typical probe is shown in Figure 3. Calibration data were 
obtained for pitch and yaw angles ranging from —20° to +20° 
in 2° increments. The Mach number of the calibration airstream 
was approximately 0.4. The independent variables of Figure 3 
are the coefficients Cp v and Cp w • Plotted together are contour 
levels of the unit vector components v and w. Solid lines are the 
Taylor’s series approximations, broken lines are calibration data. 
The lines of constant v appear nominally vertical and lines of 
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Figure 5 Cp b calibration results 
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Restrictions For Large Flow Angles 

As flow angles increase eventually this and similar algorithms 
become unreliable. There are several mathematical reasons that 
will cause the algorithm to fail. The first is when the relationship 
between the flow direction variables v and w and the flow direction 
coefficients Cp v and Cp w is no longer one-to-one. This relation- 
ship is truly one-to-one only if we restrict the set of flow direction 
variables and corresponding flow direction coefficients to within 
a neighborhood of Cp v = 0 and Cp w = 0 that satisfies Equation 
24. Where Equation 24 is not satisfied can be estimated if we 
restrict the pressure coefficients such that Cp v = Cp w = Cp m ax 
which will lead to two degenerate fourth order equations. Equa- 
tions 25 and 26. The roots of Equations 25 and 26 are determined 
using the quadratic equation. 


' Qv 



dCp* 



(24) 


(5c 4 4 C5 — 3C6)Cpt»ax + (3C2 — C^Op^ax + Cl = 0 (25) 


5(C 4 4 Cs + Ce)Cp^nax + 3(c2 + C 3 ) Cp^ciX + Ci = 0 (26) 


A second reason for the algorithm to fail is when the relation- 
ship between probe pressures p 5 and p avg and total and static 
pressure p 0 and p is no longer one-to-one. This restricts the set 
of flow direction variables and corresponding flow direction co- 
efficients to a domain containing Cp v = 0 and Cp w = 0 that 
satisfies Equation 27. Where Equation 27 is not satisfied is es- 
timated by the same restriction, Cp v — Cp w = Cpmax , which 
leads to another degenerate fourth order equation. Equation 28. 


det 


apo 

QPavq dVavq 
&P0 &p 


= ( Cp 5 - Cpav 9 f * o (27) 


(2(62 - a 3 ) + (63 - a A ))Cp* max + (2g) 

2(6l — 02 )Cpmax + 1 — Ol = 0 

All roots of Equations 25, 26 and 28 are examined and the 
smallest positive root is substituted into Equation 29. With the 
flow direction variables restricted to a domain satisfying Equation 
29 both the relationship between the flow direction variables and 
the flow direction coefficients is one-to-one and the relationship 
between probe pressures p b and p avg and total and static pressure 
is one-to-one and the algorithm will yield meaningful results. For 
the five-hole probe calibration results presented in Figures 3-5, the 
restricted domain is approximately v 2 4- w 2 < 0.25 or flow angles 
approximately less than 30°. The size of the restricted domain for 
this example is being limited by Equation 24. 

v 2 4* w 2 <2(c 1 Cp max4‘ 

(c 2 + C3)Cp 3 max + (29) 

( c 4 + C 5 -f Cs)Cp b max ) 2 


UNCERTAINTY ANALYSIS 

Because of the new algorithm’s simplicity and accuracy, it is 
well-suited to an analytical treatment (as opposed to a “jitter” 
or numerical technique) of the propagation of uncertainty in five- 
hole probe measurement The objective of the uncertainty analysis 
is to quantify uncertainty of five-hole probe measurement results 
(e.g., total pressure, static pressure, and flow direction) and de- 
termine the dependence of the measurement uncertainty on the 
uncertainty of all underlying experimental and calibration mea- 
surements. These results will provide a general procedure that 
other researchers may use to determine five-hole probe measure- 
ment uncertainty and provide guidance to improve five-hole probe 
measurement technique by pinpointing greatest sources of uncer- 
tainty. 

To discuss uncertainty, we have adopted the terminology of 
Kline (1985), where me as ur and describes a physical quantity that 
is measured and result describes a quantity that is obtained from 
a calculation involving measurands. Mathematically, results are 
dependent variables and measurands are independent variables. 
Equation 30 shows the dependence of the uncertainty of the result 
8R upon the uncertainty of the n measurands 8x \ In 
our case the results are total pressure, static pressure, and the two 
unit vector components v and w (collectively referred to as flow 
direction). 


<“>’=(l£ 4 '‘) !+ '' + (£“-)’ 001 

The uncertainty of the flow direction results calculated in data 
reduction ultimately depend on the uncertainty of the five values 
of probe pressure measured in the experiment (referred to as 
experimental pressure measurands) and the uncertainty of each 
measurement of the five probe pressures and flow angles acquired 
in probe calibration (referred to as calibration pressure measurands 
and calibration flow direction measurands). The dependence of 
flow direction results upon all intermediate steps of data reduction 
is illustrated in Figure 6. The darker paths represent values or 
results that depend wholly on experimental data and lighter paths 
values or results that depend wholly on calibration data. The 
uncertainty of flow direction results are determined from repeated 
use of Equation 30. For instance, 8v and 8w depend on 8Cp vt 
8Cp w , and 6 c\ } . . . ,6ce. The uncertainty 8Cp v depends on 
6p2,Sp\^8p b and 8p avg anti so on. This procedure is repeated 
until all the dependence shown in Figure 6 is accounted for. 

Consider the common situation when either a single pressure 
transducer is used to measure all pressures or several equivalent 
pressure transducers are used. Then, their individual uncertain- 
ties will be equal. That is, allowing for the possibility that the 
transducer(s) used for calibration and experiment are different, 
for calibration 8pi t k = ••• = 8p$ t k = 8p ca i for all k, and for 
experiment 8p x =••• = 8p b = 8p cxp . Likewise, if the same 
instrument is used to measure pitch and yaw angles to deter- 
mine flow direction throughout calibration then their individual 
uncertainties are equal, in which case 86k = &<t>k — 8>V- With 
these conditions the uncertainty of the flow direction result de- 
pends on the product of three measurand uncertainty parameters; 
8pex P /{po — p ), 8p C ai/(po — p) and 8(p , with their influence func- 
tions; / Pexp , f Pcal and /^, according to the relationship given by 
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Figure 6 Propagation of uncertainty in flow direction results 



Figure 7 Uncertainty in flow direction results 


Equation 31. Note that both measurand pressure uncertainties have 
been nondimensionalized with dynamic pressure. In Equation 31 
uncertainties for both flow direction results have been combined 
to give an overall measure of flow direction result uncertainty. 


(Sv) 2 + (6w) 2 



(31) 


The three measurand uncertainty influence functions in Equa- 
tion 31 are functions of flow direction. Plotted in Figure 7 are 


Experiment 



Figure 8 Propagation of error in 
total and static pressure results 

contour levels of the influence functions. The data used to vali- 
date the calibration and reduction algorithm were used to evaluate 
the influence functions. The contour plots of the influence func- 
tions are over the range of flow direction used when calibrating the 
probe, that is pitch and yaw angles varying from —20° to +20°. 
Several conclusions may be drawn from Figure 7. The value of 
both influence functions f Pcal and are zero when v = 0 and 
w = 0. This is not true for / Petp , however, which has the value 
/pci p (0,0) = 2ci/(l — ai). This is a general observation that 
is valid for all five-hole probes calibrated with this algorithm. In 
practical terms this means that for modest flow angles, most or 
all flow direction result uncertainty comes from uncertainty in the 
experimental pressure measurands. For large flow angles, near 
the boundary of the calibration domain, the values of and par- 
ticularly f Pcal grow rapidly. However, levels of f Pcal and 
generally are less than f Pcxp except for large values of both v and 
w where f Pcxp actually decreases. 

Like flow direction, the uncertainty of total and static pressure 
results calculated in data reduction depend on the uncertainty of 
the five values of probe pressure measured in the experiment, 
the uncertainty of each measurement of the five probe pressures 
acquired in probe calibration as well as each measurement of the 
calibration stream’s total and static pressure. Figure 8 shows the 
dependence of total and static pressure results on intermediate 
data reduction results. Because flow direction does not appear 
in Figure 8 the uncertainty of total and static pressure results 
do not explicitly depend on calibration flow direction measurand 
uncertainty. Values of total and static pressure result uncertainties 
are determined by repeated use of Equation 30. 

Figure 9 shows contour levels of the two measurand uncertainty 
influence functions for total pressure in Equation 32. There are 
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Figure 9 Uncertainty in total pressure results 

several similarities between the total pressure result uncertainty 
influence functions and the corresponding flow direction result 
uncertainty influence functions shown in Figure 7. First, the value 
of g PcaJ is zero when v = 0 and w = 0 while g Pcxp ( 0,0) is 
nonzero. Like flow direction result uncertainty, for small flow 
angles the total pressure result uncertainty comes almost entirely 
from uncertainty in the experimental pressure measurands. In fact, 
g P cx P (0,0) = 1 which simply means that it is impossible for total 
pressure result uncertainty to be less than experimental pressure 
measurand uncertainty. The value of g Pcal increases quickly away 
from the origin, however, values of g Pcal remain less than g Pcxp . 



Figure 10 shows contour levels of the two measurand uncer- 
tainty influence functions for static pressure result uncertainty in 
Equation 33. There are several interesting differences with the 
corresponding influence functions for flow direction result uncer- 
tainty and total pressure result uncertainty. First, the value of 
h Pca i is nonzero when v = 0 and w = 0, while h Pcxp (0, 0) = 1. 
This means that for small flow angles, the uncertainty of static 
pressure results are necessarily greater than the uncertainty of to- 
tal pressure results. This fact confirms the belief that static pres- 
sure determined by a five-hole probe is not as “accurate” as total 
pressure obtained with the same probe. 



The preceding analysis can be used to provide guidance in ways 
to reduce five-hole probe result uncertainty. In general, reducing 
result uncertainty can be achieved by either reducing the measur- 
and uncertainty parameters 6p cxp /(po — p), 6p ca i/(po — p) and 
6<p or by reducing their respective influence functions. The mea- 
surand uncertainty parameters 6p cxp /(po — p) and 6p C ai/(po — p) 
can be reduced by either decreasing the numerator by using more 
accurate transducers or by averaging repeated measurements, or 


by increasing the denominator by raising the dynamic pressure of 
the flow. However, the ability to increase dynamic pressure is 
limited since this typically requires higher velocities and as the 
sonic velocity is approached compressibility effects become sig- 
nificant The measurand uncertainty parameter 6<p can also be 
reduced by using a more accurate angle measuring device or av- 
eraging repeated measurements. 

The measurand uncertainty influence functions / Pexp , y Pexp and 
h Pcxp depend on the probe calibration coefficients as well as v 
and w. To reduce f Pcxp% g Pcxp and h Pcxp requires changes in the 
probe calibration coefficients, which are determined by the probe 
geometry. To reduce f Pcxp , g Pcxp and h Pcxp requires making 
physical modifications to the probe. For flow direction uncertainty, 
fpexpi 0>0) = 2ci/(l — ai). Using the approximations stated 
earlier relating probe calibration coefficients to probe geometry 
yields the approximation f Pcxp ( 0,0) « 0.5/ sin Therefore, 
flow direction result uncertainty can be improved by increasing 
the angle of the openings of the four outer tubes. However, 
intuition says that this will also restrict the range of flow angles for 
which the probe can be used. For total and static pressure result 
uncertainty y Pexp (0,0) = /i Pexp (0, 0) = 1. This observation is 
valid for all five-hole probes, so it is not possible to lower this 
value by making modifications to the probe. Modifying the probe 
will not improve total or static pressure result uncertainty for small 
flow angles. 

The measurand uncertainty influence functions / Pcal , y Pcai , 
h Pcal and depend on calibration data as well as v and w. The 
probe calibration coefficients themselves to not explicitly appear 
in f Peal , g Pcal , h Pcal , and / v , however modifications to the probe 
geometry will effect them by changing the calibration measure- 
ments. The influence functions f Pcal , g Pcan /i Pcal , and may 
be reduced by modifying the probe calibration procedure. One 
way of accomplishing this is by acquiring more calibration data. 
For n total calibration measurements used the influence functions 
Lr 9pc i. Kc * ^ ft ^ approximately 0( 1/y/n) so the 
overall levels of f Pcal , g Pcal , h Pcal and may be approximately 
halved by acquiring four times as much calibration data. To ver- 
ify this the uncertainty values presented earlier were compared to 
uncertainty values calculated using a subset of the original cal- 
ibration data containing one-fourth of the number of calibration 
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Figure 11 The Internal Fluid Mechanics Facility (IFMF) (Dimensions in feet) 


data. This subset was obtained by skipping alternating rows of 
pitch and yaw data so that the distribution of pitch and yaw an- 
gles for the calibration data subset was equivalent to the original 
calibration data set. This analysis showed that the overall values 
of the influence functions f Peal , g Pcal , h Peal and were nomi- 
nally greater by a factor of two for the smaller calibration data set. 
The influence functions f Pcxp > g Pcxp and h Pcxp were unchanged, 
as was expected 

Changing the distribution of pitch and yaw values used for 
calibration will also change f Pcal , g Pc ml , h Pcal * and / v . In 
general, these influence functions are reduced when the pitch 
and yaw values used for calibration are concentrated at higher 
flow angles. This was confirmed by using two subsets of the 
original calibration data. Both subsets contained the same amount 
of data, but one subset used calibration data acquired at large 
flow angles while the other used data acquired at small flow 
angles. The resulting values of / PcaI , g Pcat > » and /<, were 
compared and showed interesting results. For total pressure result 
uncertainty, the value of g Peal at large flow angles from the large 
flow angle calibration data was significantly less (by a factor of 
four) than the corresponding small flow angle calibration data set, 
while g Pcal (0, 0) = 0 for both calibration data subsets (as stated 
earlier). The same observations hold for f Pcal and / v . For flow 
direction result uncertainty and total pressure result uncertainty the 
conclusion is unambiguous. Given a fixed amount of calibration 
data, flow direction and total pressure result uncertainty is reduced 
by using calibration data acquired at larger flow angles. For static 
pressure result uncertainty the comparison of h Pcal was more 
complicated. Values of h Pcal at large flow angles were less for 
the large flow angle calibration data set, as was expected, but 
h Pcal (0, 0) was less for the small flow angle calibration data set. 
This suggests that static pressure result uncertainty is improved in 
the region where the calibration measurements are concentrated. 

The major results of the uncertainty analysis are summarized 
below: 

1. The uncertainty of five-hole probe results (e.g., total pres- 
sure, static pressure and flow direction) are shown to depend 
on the product of three measurand uncertainty parameters 


and their corresponding influence functions. The measur- 
and uncertainty parameters quantify the uncertainty of the 
underlying calibration and experimental measurands. 

2. For modest flow angles, most or all flow direction result 
uncertainty and total pressure result uncertainty comes from 
uncertainty in the experimental pressure measurands. 

3. Total and static pressure result uncertainty can not be less 
than experimental pressure measurand uncertainty. 

4. For modest flow angles, static pressure result uncertainty is 
necessarily greater than total pressure result uncertainty. 

5. Result uncertainty may be reduced by reducing either the 
measurand uncertainty parameters or their influence func- 
tions. The measurand uncertainty parameters may be re- 
duced by either using instruments with less measurement 
uncertainty or by using repeated measurements. 

6. To reduce the three measurand uncertainty influence func- 
tions corresponding to experimental pressure measurand un- 
certainty requires modifying the probe geometry. For modest 
flow angles, flow direction result uncertainty can be reduced 
by increasing the angle of the openings of the four outer 
tubes. For modest flow angles, the total and static pressure 
result uncertainty coming from experimental pressure mea- 
surand uncertainty can not be reduced. 

7. To reduce the four measurand uncertainty influence functions 
corresponding to calibration pressure and flow direction mea- 
surand uncertainty requires modifying the calibration proce- 
dure. The influence functions can be reduced by approxi- 
mately one-half by increasing the amount of calibration data 
acquired by a factor of four. 

8. In general, flow direction and total pressure result uncertainty 
may be reduced by acquiring calibration data at larger flow 
angles. However, this is not generally true for static pressure 
result uncertainty. 

THE FIVE-HOLE PROBE RAKE 
Motivation 

The five-hole probe rake is built for use in NASA Lewis* 
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IFMF. This facility is schematically illustrated in Figure 11. The 
test section Mach number range is between 0.0 and 0.8 with 
corresponding mass flow rates between 0 and about 16 lbs/sec. 
Details concerning the operations and capabilities of the IFMF 
may be found in Potto et al. (1991). The IFMF is currently 
being used to test various inlet diffuser designs. These diffuser 
geometries possess high degrees of stream wise curvature and large 
•changes in cross-sectional area providing a catalyst for strong 
secondary flows and attendant boundary layer flow separations. In 
the past, the resolution of these three-dimensional velocity fields 
was accomplished with a single five-hole probe traversed in the 
cross plane of the test section exit duct. The time required for a 
complete survey was considerable. As a means of saving time, 
the five-hole probe technique was extended to a rake arrangement; 
this idea prompted by the common use of Pitot probe rakes for 
measurements of total pressure. 

Geometry, Construction, and Calibration 

The rake probe is illustrated schematically in Figure 12. Ten 
individual five-hole probe tips are mounted one inch apart on 
the cylindrical stem spanning the tunnel diameter. Each tip is 
constructed from five 0.020” outer diameter stainless steel tubes 
silver brazed in the pattern indicated in Figure 12. After brazing, 
each tip face is ground to a 45° conical shape using a lathe and 
a fine grinding stone. A nominal tip diameter of 1/16” results 
from the brazing and grinding process. Each tip is individually 
calibrated following the previously outlined procedure. Probe 
symmetry is gauged by inspecting the calibration curves like those 
shown in Figures 3 through 5. Probe tips not meeting the desired 
level of symmetry are reground and recalibrated, or discarded. 



Once mounted in the stem, the relative misalignment in pitch and 
yaw between consecutive tips is measured by traversing the full 
length of the rake through the calibration facility’s open jet and 
recording the response of each tip. 

Application and Results 

Figure 13 illustrates the position of the rake probe inside the 
exit duct One end of the rake probe is secured to a linear motion 
actuator which provides probe traversing in the radial coordinate, 
r. The other end of the probe also passes through the duct wall, 
and through a device similar to a linear bearing. This arrangement 
greatly reduces probe vibration and deflection when compared to 
the method of cantilevered support used with conventional probes. 
This can be a significant benefit when flow velocities are large or 
the distance to be spanned in a survey is large. Traversing in the 
circumferential coordinate, (, is accomplished by manually turning 
the duct with the rake probe and actuator attached. Surgical grade 
TVgon tubing (0.017” ID) is routed from the aft end of each probe 
tip, out of the exit duct, and to a series of Electro Scanning 
Pressure (ESP) transducer modules located in the facility. Probe 
pressures are recorded by NASA Lewis’ ESCORT data acquisition 
system. The collected pressure data are transferred from storage 
on the ESCORT Scientific VAX cluster to Sun workstations where 
the reductions to velocities are performed. 

Figure 14 illustrates transverse velocity components as deter- 
mined by the rake probe at the exit plane of an S-duct, which is 
the current inlet diffuser model undergoing testing at the EFMF. 
The grid resolution is Ar = 0.25*, AC = 10°. The local Mach 
number is approximately 0.4 and the total conditions are approxi- 
mately atmospheric. Figure 15 is the corresponding result obtained 
with a single five-hole probe (Wellborn et al ., 1992). The only 
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Figure 14 Transverse velocity components at the diffusing 
S-duct exit plane as measured by the rake probe 


significant difference between Figures 14 and 15 is the time re- 
quired to collect the data. The data illustrated in Figure 14 was 
collected in less than 5% of the time required for Figure 15. Use 
of the rake allows a substantial reduction in time and effort for 
the researchers at NASA Lewis’ IFMF. 



Figure 15 Transverse velocity components at the 
diffusing S-duct exit plane as measured by a 
conventional five-hole probe (yaw-nulling method) 


SUMMARY 

The large data sets associated with the use of a rake of five- 
hole probes provides the motivation for a novel algorithm of five- 
hole probe calibration and data reduction. In this algorithm four 
pressure coefficients are defined and related to flow angles through 
Taylor’s series. Probe tip symmetries are used to reduced the 
number of series coefficients from sixty to thirteen. Laboratory 
calibrating and testing of individual five-hole probe tips verifies 
the validity and accuracy of the algorithm. In data reduction, 
the algorithm allows a much faster conversion of pressure to 
velocities, making the rake arrangement of five-hole tips a practical 
expedient to acquiring three-component velocity data sets in duct 
flow research. 

An analytical treatment of the propagation of uncertainty in 
five-hole probe measurement was conducted. The uncertainty 
analysis determined the dependence of the result uncertainty on 
the uncertainty of all underlying experimental and calibration 
measurands and revealed details about the relationships between 
probe geometry, calibration, data reduction and result uncertainty. 
This study outlines a general procedure that other researchers may 
use to determine five-hole probe result uncertainty and provides 
guidance to improve measurement technique. 

The ability to measure total and static pressure (hence velocity 
magnitude) and flow direction with a five-hole probe is well 
established. Recently we have constructed a rake of five-hole 
probe tips that operates in a non-nulling mode. Current non- 
nulling five-hole probe calibration and data reduction procedures 
typically involve complex interpolate algorithms. Development 
of the new algorithm makes the rake probe practical to use. The 
new algorithm uses two components of the unit vector in the flow 
direction instead of pitch and yaw angles as variables as well as 
probe tip symmetries to greatly simplifying both calibration and 
data reduction. In data reduction, total pressure, static pressure, 
and flow direction are calculated directly, without tables, spline 
network interpolation, or additional approximations. The result 
is a much faster reduction to velocities, making a rake probe 
arrangement of five-hole tips a practical expedient in duct flow 
research. 
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APPENDIX CALIBRATION AND DATA REDUCTION 

The coordinates of the ith probe opening (measured relative 
to the probe coordinate system) is x< = (z,, y,, a,-). Because of 
probe zz-plane symmetry the following transformations are valid 

T »( x i) = x i 
T,(x 2 ) = x 4 
T s (x 3 ) = x 3 
Tj,(x 4 ) = X 2 
T y( x s) = x 3 

where T y : (x, y, z) — ► (x, —y, z ). 1 For xy-plane symmetry the 
following transformations are valid 

T z (xi) = x 3 
T*(x 2 ) = x 2 
T z (x 3 ) = xi 
T z (x 4 ) = x 4 
T z (x s ) = x 5 

where T* : (x, y, z) — * (x, y, —z). Likewise, for 90 ° rotationsi 
symmetry about the x-axis the following transformations are valid 

T9o(Xl) = x 4 
T 90 (X 2 ) = X! 

T90 (x 3 ) = x 2 
T 90 (x 4 ) = x 3 
T 90 (x 5 ) = x 5 

where T 90 : (z ,y,z) — ( x,z,-y ). 

Consider a five-hole probe in a steady, uniform air flow. The 
boundary conditions that define the flow field are the free stream 
conditions at a distance well upstream of the probe. They are 
the free stream velocity Voo = {Uoo, Voo, Woo) (measured rel- 
ative to the probe coordinate system) and the free stream pres- 
sure poo . For given free stream conditions there is a functional 
relationship between spatial coordinates and velocity and pres- 
surre, V : (x,y, z) -*• ( U,V,W ), P : ( x,y,z ) —■ p and 
P 0 : (z,p, z) — • - po. The pressure sensed by the ith probe open- 
ing is pi = p(x, ) 

Because of zz-plane symmetry and the invariance of the equa- 
tions of motion T s oVo T y 2 and P o T y are the velocity and 
pressure fields defined by the free stream conditions T y o Voo 
and poo • This requires 

Pi = p(xi) = P o T y (xi) 

P2 = p(x 2 ) = P o T y (x 4 ) 

Pz = p(x 3 ) = P o T y (x 3 ) 
p 4 = p(x 4 ) = P o T y (x 2 ) 

Ps = p(x 5 ) = P o T y (x 5 ) 


1 The function T y is sometimes called a flip. The domain of T y 
is not necessarily physical space. For instance T y : (ti, v, w) — ► 
( w , —v, u>)« 

2 The notation o used here is the composition of two functions. 
For instance, for functions f{x) and g{x) then / o g = f{g{x)). 


The expression PoT y (x,) is interpreted as p t for the free stream 
conditions T y o Voo and po©. For instance p 2 for free stream 
conditions Voo and poo equals p 4 for free stream conditions 
T y o Voo and poo- Thinking of the probe pressures as functions 
of free stream conditions then 

Pi {Uoo, Voo , W©© , Poo) = Pi {Uoo , — Vx> , Woo j Poo ) 

P 2 (f^oo , Voo > W©©, Poo ) = Pi{U oo, Vqo , Wqo , Poo ) 

Pz(Uoo, Voo , Woo 1 Poo) = P 3 {Uoo, Voo , Woo , Poo ) 

Pi{Uoo, Voo, Woo,Poo) = P 2 (f^OO j Voo , W©© , Poo ) 

Pb {U OO, Voo, Woo, Poo) = Pb {Uoo, Voo, Woo, Poo) 

For zp-plane symmetry T z o V o T z and ?oT : are the 
velocity and pressure fields defined by the free stream conditions 
T z o Voo and poo- Therefore 

Pi = p(xi) = P o T z (x 3 ) 

P2 = p(x 2 ) = P o T z (x 2 ) 

Pz = p(x 3 ) = P o T*(xi) 

P4 = p(x 4 ) = P o T z (x 4 ) 

Pb = p(x 5 ) = P o T z (x 5 ) 

and 

Pl(UoO,Voo,Woo,Poo) = Pz(U 00, Vqo, — Woo, Poo) 

P 2 (Uoo , Voo , Woo , Poo ) — P 2 {Uoo, Vx), W©© , Poo) 

PZ (Poo , Vqo » Woo , Poo ) = Pi (Uoo , Voo , — Woo, Poo) 

Pi {Poo, Vqo , Woo , Poo ) — Pi {Uoo, Voo , Woo, Poo) 

Pb{U 00 , Voo , Wqo , Poo ) = Pb{U 00 , Vqo , — Woo, Poo) 

For 90 ° rotational symmetry T90 o V o T90 and P0T90 are the 
velocity and pressure fields defined by the free stream conditions 
T90 o Vqo and poo- Therefore 

pi = p{xi) = P o T 90 (x 2 ) 

P2 = p(x 2 ) = P o T90 (x 3 ) 
p 3 = p(x 3 ) = P o T 90 (x 4 ) 
p 4 = p(x 4 ) = P o T 90 (xi) 

Pb = p(x 5 ) = P O T 9 o(x 5 ) 

and 

Pi {Uoo, Vqo , Woo , Poo) = P2 {Uoo , W 00 , — V»,Poo) 
MUoo^W^Poo) = p^Uoo.Wco , Vqo, Poo) 

PZ {Uoo, Vqo, Wqo, Poo) = P 4 {Uoo , Wqo , Vx>»Poo) 

Pi {Uoo, Vqo , Wqo , Poo ) = Pi {Uoo, Wqo , V©©, Poo) 

Pb{U 00, Voo , Wo© , Poo ) = Pb{U oo, Wqo , Vqo, Poo) 

To develop our calibration and data reduction procedure, five 
pressure coefficients, one for each probe opening, are defined by 
Equation 1 . These coefficients are not used for calibration and data 
reduction, but they help analyze concepts relating to probe symme- 
try that are essential to our calibration and data reduction scheme. 
Since p t may be considered a function of U Q o , Vqo , Woo , Poo then 
so is the coefficient Cpi a function of these four quantities. Like 
the pressure coefficients described previously the pressure coeffi- 
cients given by Equation 1 are insensitive to a range of Reynolds 
and subsonic Mach numbers and depend on flow direction only. 
Therefore, Cp, is a function of the unit vectors Uqo, Vqo, tt>oo- 
Since this is a unit vector, only two components are independent, 
so with the restriction Uoo > 0 and dropping the 00 subscript for 
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the sake of simplicity we say Cp, is a function of v, to, as is 
shown in Equation 1. 


Cpi = 


P\{U oo •> Voo ) OO ) Poo) Po 


PO.oo — Poo 

= Cpi(v,w) 


Cp v (v,w) = -Cp v (-t>, to) 
Cp v (v,to) = Cp v (v, — to) 
Cp w (v, to) = Cp w (—v, to) 
Cp w (v y w) = — Cpxu (v , —to) 

(1) Cp v (v,w) = Cp w (w, -v) 

Proofs for these rules are stated below 


We may now state the probe symmetry rules in terms of the 
pressure coefficients Cp,. For zz-plane symmetry 

Cp\(v,w) = Cpi(-v,w) 

Cp 2 (v, to) = Cp\(—v,w) 

Cp 3 (v,to) = Cp 3 (-v,to) 

Cp 4 (v,to) = Cp 2 (-v,to) 

Cp 5 (i>,to) = Cp 5 (-v,to) 

For xy-plane symmetry 

Cpi(v,to) - Cp 3 (v,-to) 

Cp 2 (v,to) = Cp 2 (v, —to) 

Cp 3 (t>, to) = Cp\ (v, —to) 

Cp 4 (v,to) = Cp 4 (v,-to) 

Cps(v, to) = Cp$(v, —to) 

For 90° rotational symmetry about thex-axis 

Cp\(v,w) = Cp 2 (w,-v) 

Cp 2 (v,to) = Cp 3 (to,-v) 

Cp 3 (t?, to) - Cp\ (to, —v) 

Cp 4 (v, to) = Cpi(w,-v) 

Cp 3 (v,w) = Cps(w, —v) 

Proofs of the rules stated above are very straightforward. For 
instance the proof that Cpi(t>, to) = Cpi(— v,to) is 

r>, ( \ Pi (^oo, Vqo, ^<x> , Poo) —poo 

Cpi(v,w) = 

P0,oo Poo 

PljUpO) Voo, OO ) Poo ) Poo 

p0,oo — Poo 

= Cpi(-v,w) 

The remaining proofs are nearly identical and are not repeated. 

The pressure coefficients used in calibration and data reduction 
and given by Equations 1, 3 and 4 can also be defined in terms of 
the five pressure coefficients, as shown below. 


Cpavg( V, W) = ^(Cpi (V, w) + Cp 2 (v, w) 

+ Cp3(v, w) + Cpt(v , «;)) 

= t(<?Pi (-»> «0 + Cp4 (-' W. w) 
4 

+ Cp 3 (-v,to) 4-Cp 2 (-t;,to)) 
— Cp<»vy ( — v, to) 


Cpavg(?, W) = ^(Cpi(v>w) + Cp 2 (v y U>) 

-F Cp 3 (v, to) 4- Cp 4 (v, to)) 

= ^(Cp3(v, -to) 4- Cp 2 (v, -to) 
4- Cpi (v, -to) + Cp 4 (v, -to)) 

= Cp<iv g(v , — to) 


Cpavg(v,w) = i(Cpi(v,U;) + Cp 2 (v,w) 

4- Cp 3 (o, to) 4- Cp 4 (v, to)) 

= 7(^2 (to, -t>) 4- Cp 3 (to,-v) 
4 

+ Cp 4 (to, -v) 4- Cpi (to, -v)) 

= Cpavy(tO,-v) 


Cp v (v, to) 


Cp 2 (t;, to) - Cp 4 (u, to) 
Cp 5 (t>, to) - Cp 

av y(v, to) 

Cp 4 (-o, to) - Cp 2 (— o, to) 
Cp 5 (-t>, to) - Cp OV y(-V, to) 
Cp2(-t;, to) - Cp 4 (— o, to) 
^Ps(-V, to) - C7p av9 (-V, to) 
-Cp v (— 0, to) 


Cp<xvg = T^Pl 4" Cp 2 4* Cp 3 + Cp 4 ) 
4 


Cp v — 


Cp w = 


Cp 2 — Cp 4 

Cps — Cpavg 
Cjp 3 — Cpi 


Cps — Cpavg 

The probe symmetry rules for the coefficients Cp av9 . Cp v and 
Cpw are 

Cp av9 (v,t0) = Cpav 9 (-V,to) 

Cp at> 9 (t;, to) = Cpavg(v,-w) 

Cpavg(v,w) = Cpat/ 9 (t 0 , -v) 


Cp v (u, to) 


Cp 2 (u,to) - Cp 4 (t;, to) 
Cp 5 (v, to) - Cp aV9 (v, to) 
Cp 2 (u, -to) - Cp 4 (u, —to) 
Cp 5 (v, -to) - Cpav 9 (v, -to) 
£Pv(t>, -to) 


Cpu,(v, to) 


Cp 3 (y, to) - Cpi(v, to) 
Cp 5 (v, to) - Cpavg(v , to) 
Cp 3 (-U,t0) ~ Cpi(-U, to) 
Cp 5 (-V,to) - Cpav 9 (-V,to) 
Cpu;(-V>to) 
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Cp w (y , w) 


Cpz(v,w) - Cpi(v,w) 
Cp$(v,w) - Cpavg(v,w) 

Cpi ( V , -ty) - Cp 3 ( V, -w) 
Cp 5 (t>, -tl>) - Cpav<7(v, -w) 
Cp 3 (t;, -w) - Cpi (v, -ty) 

Cp 5 (v, - w) - Cpavg(y , -ty) 

— Cp w (»,— to) 


Cp,,(v,ty) 


Cp2(ty,u;) ~ Cp 4 (v, w) 

C^ 5 (v, ty) - Cp av g(V t W) 

Cp 3 (m,~ t>) - Cpi(tD,-tl) 
Cps(w, -v) - Cpavg{w, -v) 
Cp w (w,-v) 


The symmetry rules for v and w as function of Cp v and Cp w 
are stated and proven below 


v(Cpv)Cpw) — v( CpviCpw') 
v{Cpv j Cpw') — v(C7pt/ , “ Cpw ) 
w(Cp v ,Cpu,) = ty(-Cp„,Cpu,) 
w{^C/Pxj f Cpxv') ~ w(Cp v , Cpu;) 
»(CPv, Cpw) = ty(Cpu;, — Cp V ) 


v(Cp v ,Cpu,) = v(Cp v (t;,ty),Cpu,(t;,iy)) 

= t>(Cp„(t>, -ty), Cp,x,(v, -ty)) 

= v(Cp v (v , ty), -Cp w (t>, ty)) 

= t>(Cp v , -Cpu,) 

The proof that w(Cp v ,Cp w ) — w(—Cp Vi Cp w ) is identical and 
is not repeated. 

v(Cp v ,Cp w ) = v(Cp v (v 1 w),Cp w (v,w)) 

= -v(Cp v (-t;, ty), Cp w (-v, ty)) 

= -v(-Cp v (t;, ty), Cp w (v, ty)) 

— v( Cp v1 Cp w ) 

The proof that w(Cp v ,Cp w ) = — ty(Cp„, — Cp w ) is identical 
and is not repeated. 

v(Cp v (v , ty), Cp w (y , ty)) = w(Cp v (w, v), Cp w (w , v)) 

v(Cp v ,Cp w ) = ty(Cpu,(v, -ty), Cp„(t;, -ty)) 

= ty(-Cpu,(v, w),Cp v (v,w)) 

— - ty( Cp w , Cpv ) 

= w(Cp w , Cp v ) 

The symmetry rules for Cp a vg and Cps as functions of Cp v 
and Cp w are stated and proven below. 

Cpavg^Cpv , (sPw) — Cpat/g( C'PviCpw') 

Cpav gi^Cpv i (-'Pw') = Cpavg(Cp v , Cpu,) 

Cpavg(Cp v , Cp w ) = Cp a vg{Cp w , — Cp„) 

C/ps(Cpv i Cpw^ = Cps( CPvj^Pti;) 

Cps {Cpv , Cp w ) = Cps(Cpvi Cp^) 

Cp5 (C'pv, C/Pu,) = Cps^Cpwi Cp v ) 


Cp a vg{CPv , Cpw') — Cpavg{Cpv (tJ, ty), Cp w (t 7 , ty)) 

= Cpavg(Cp v (—v, ty), Cp w (— 1>, ty)) 

= Cpavg(-Cp v (v , ty), Cpu,(v, ty)) 

= C , p OV 5( — Cpvi Cpty) 

The proofs that Cpovj(^Pv> Cpw) = (-^Pavg^^Pv , ~~Cp u;)» 
Cps{Cpv , Cpw ) = Cps( ^PvjCpty) and Cp5(Cp v , C7pty) — 
Cp 5 (Cp v , — Cpu,) are identical and are not repeated. 


Cpav g(Cpv , Cpty) — Cp a t/g(Cpv (v, ty), C^Pty (t7, ty)) 

= Cp aV 9 (Cpu,(ty, -v), Cp v (ty, -t>)) 
= Cp a vg(Cpw i Cpv ) 

— Cpavg(Cpwt Cpv) 


The proof that Cps(Cp v , Cp^) = Cps(Cp w , —Cp v ) is identical 
and is not repeated. 

The consequence of the symmetry rules on the Taylor’s series 
for Cpavg and Cps are stated and proven below 

a *j — t = 1,3,5, ... 

a *j = j = 1,3,5,... 

a ij = a ji 

bij = 0, t = 1,3,5, ... 
bij = 0, J = 1,3,5,... 

bij — bji 


Cpav g[CpV) Cpw') — ^00 

+ QloCpv + CL 0 lCp w 
+ CL 20 Cpi -f 0 L\\Cp v Cp w + CL 02 Cp 2 w 
+ a>3oCpl + a2\Cp 2 v Cp w + CL^CpvCpij 

+ <*03 Cp* w 

+ <*4oC'Pv + 0L3\Cp Z v Cp w 4- Q’22Cp 2 Cp^jj 
+ <*i 3Cp v Cp z w -f ao4 Cp K w 
+ ... 


Cpavg( Cp v ,Cp w ) — aoo 

— CLioCp v + aoiC'pu, 

+ d2oCpl — anCpt/Cpu, 4- <*02 Cp^, 

— <*30 CpJ 4* <*21 Cp 2 Cp w — 0.\2Cp v Cp w 
4 aosCp^ 

+ <*4oC^t ~ 4* <*22 Cp 2 Cp 2 w 

— <* 13 CpvCpw 4- <*04 Cpi, 

+ ... 
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Cp<iv g{Cpv > Cp w ) Cpavg( Cp»i Cp w ) = 

2dioCp v 

4~ 2ct 1 1 C7 pvCpw 
+ 2a3oC’j?v 4- 2ai2C'p v C , p^ ; 

4* 2031 ( 7 ^^*?^ 4~ 2ai3 Cp v Cpl 

+ ... 

= 0 

Since Cpav g(C p%> f Cp%v) ^Povj( Cpt), Cpiu) — 0 for all Cpv 
and Cp w then all the coefficents a,j where i = 1,3,5,... are 
equal to zero. The proofs that a tJ =0, j = 1, 3, 5, . . ., bij = 
0, t = 1, 3, 5, . . . and bij = 0, j = 1, 3, 5, ... are identical and 
are not repeated. 

Cpavg(Cp v ,Cp w ) = <*00 

4- aioC'pv + <*01 Cp w 

4- <*20 Cp 2 4- o.\\Cp v Cp w 4- ao2Cp^ 

4- <*30 4- <*21 Cp 2 Cp w 4- <*12 Cp v Cpl 

4- &ozCp w 

4” CL4oCp* 4- dZ\Cp Z v Cp w 4~ 0'22Cp 2 v Cp\ ) 

4- CL\zCp v Cpl 4- <*04 Cpl 
4-... 

Cpavg{Cpw > Cpv ) = <*00 

4- &ioCp w 4“ <*01 Cp v 

4- OL2oCp 2 w 4- 0L\\Cp w Cp v 4- <*02 Cp 2 

4- <*30 Cpl 4- <*21 Cp 2 w Cp v 4- <*12 CpxjjCp^ 

4- <*03 Cpl 

4- <*40 Cpl + a>Z\Cp Z w Cp v 4- 0>22Cp 2 w Cp 2 
4* &izCp w Cpl 4- Qo+Cp* 

4-... 

Cpavg(Cp v , Cpxv ) — Cpavg{Cp w , C7p v ) 

= 4-(<*io — 0 L 0 l)(Cp v — Cp w ) 

4- (<*20 — <*02 )(Cpl — Cpl) 

4- (030 — <*o z){Cpl — Cpl ) 

4- (<*21 — <*12) {Cp 2 Cp w — Cp v Cpl) 

4* (<*40 — &o\){Cp\ — Cpl) 

4- (<*31 - aiz){CplCpu, - Cp v Cp 3 w ) 

+ ... 

Since Cpav g(Cpv , Cp%v) Cpav g(Cpw , Cpv) — — 0 for all Cp v 

and Cp w then this requires aij = aji . The proof that bij = bji 
is identical and is not repeated. 

The consequence of the symmetry rules on the Taylor’s series 
for v and w are stated and proven below 

Cij — 0, * = 0,2,4, ... 


Cij — 0, j = 1,3, 5,... 
dij — 0, x — 1,3,5,... 
dij = 0, j = 0,2,4, ... 

c *i = dji 

v(Cp Vy Cp w ) — Coo 

4- cioCp v 4- coiCp w 
4- C 2 oCp 2 4- c\\Cp v Cp w 4- coiCpl 
4- czoCpl 4- C2\Cp 2 Cp w 4- c\2Cp v Cpl 
4- cozCpl 

4- C4oCpt 4- cz\Cp z v Cp w 4- C 22 Cp\Cpl 
4- cizCpvCpl 4- co\Cpl 
+ ... 

v( — Cp v ,Cp w ) — coo 

— cioCp v 4- co\Cp w 
4- C 2 oCp 2 — c\\Cp v Cp w 4- co 2 Cpl 
— C3oCpv 4- C2\Cp 2 Cp w — c\2Cp v Cpl 
4- cozCpl 

4- c\oCp A v — cziCplCp w 4- C22Cp 2 Cpl 
— cizCp v Cpl 4- co+Cpl 
4-... 

v{Cp v , Cpuf) 4" v( — Cp v , Cpxv) = 

2 coo 

4- 2coiC7pu; 

4- 2c2oCpl 4- 2co2C'p^ 

4- 2 c 2 \Cp 2 Cp w 4- 2cazCpl 
4- 2c4oCpt 4- 2c22Cp 2 Cpl 4- 2co4Cpi, 

4-... 

= 0 

Since v(Cpt;, Cpm) 4~ v(~ ■ 'Cp^^ Cpxv) — 0 for all Cp^ and C p^ 
then all the coefficents dj where i = 0,2,4,... are equal to zero. 
The proof that dij = 0, j = 0 , 2 , 4 , ... is identical and is not 
repeated. 

v{Cp v , Cp w ) — c oo 

4- cioCp,, 4- co\Cp w 
4- C 2 oCp 2 4- c\\Cp v Cp w 4- C 02 Cpl 
4- czoCpl 4- C2\Cp 2 Cp w 4- c\2Cp v Cpl 
4- cozCpl 

4- c\qCp\ 4- cz\Cp z v Cp w 4- C 22 Cp^Cpl 
4- cizCp v Cp 3 w 4“ caCpl 
+ ... 
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v{Cp v , Cpxv') — coo 

4 CloCpv — CoiCp w 
4 C 20 Cp v c\\Cp v Cp w 4 C 02 CP XV 
4- C 3 oCpv — C 2 \Cp 2 v Cp w 4 c 12 Cp v Cp 2 w 

— cozCp z w 

4 C\oCp\ — Cz\Cp Z v Cp w 4 C22Cp 2 v Cp 2 w 

— c\zCp v Cp z w 4- c<nCp* w 

+ ... 

v(C pv > Cpw') v{Cp v , Cpxv} ~ 

4- 2coi Cp w 
4* 2cn Cp v Cp XV 
4 2c2iCp 2 Cp w 4- f 2’CozCp z w 
4* 2c3i Cp z Cp w 4* 'Zc\zCp v Cp 2 w + 

+ ... 

= 0 

Since v(Cp v , Cp w ) — v(Cp Vl —Cp w ) = 0 for all Cp v and Cp w 
then all the coefficents cij where j = 1, 3 , 5 , ... are equal to zero. 
The proof that dij =0, j = 0, 2, 4, ... is identical and is not 
repeated. 

v(Cp v , Cp w ) — coo 

4* cioCp v + coiCpxv 

4" C2o^Pv 4- c\\Cp v Cp w 4- co2C7p^ 

4- C3oCpl 4- C2\Cp 2 Cpxv 4- C\ 2 Cp v Cp 2 w 
4- cozCp z w 

4- c+oCpl 4- cz\Cp z v Cp w 4- C 22 Cp 2 Cp 2 w 
4- c\zCp v Cp z w 4 co\Cp* w 
4-... 


w(Cp w ,Cp v ) = d 0 o 

4- dioCp XV 4 doiCp v 
4 d2oCp 2 w 4 d\\Cp w Cp v 4 do2Cp 2 
4 dzoCp z w 4 d2\Cp 2 w Cp v 4 d\2Cp w Cp 2 v 
4 dozCp z v 

4 d\oCp A w 4 dz\Cp z w Cp v 4 d 2 2 Cp 2 w Cp 2 
4 dizCpwCpl 4 do+Cpl 
4 ... 


v{Cpv , Opw) — w(Cpxvi Cp v ) 

= 4(cio — doi)(Cp v — Cp w ) 

4 (C20 — do 2 )(Cpl — Cplj ) 

4 (C30 — doz)(Cpl — Cp^,) 

4 (C21 — d\ 2 ){Cp 2 Cp w — CpvCplj) 

4 (C40 — do 4 )(Cpi —Cp\v) 

4 (031 — d\z){Cp\Cp w — CpvCplj ) 

4 ... 

Since v(Cp v , Cp w ) — w(Cp Wl Cp v ) = 0 for all Cp v and Cp w 
then this requires c tJ = dji. 


18 


REPORT DOCUMENTATION PAGE 

Form Approved 
OMBNo. 0704-0188 

Public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, 
gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this 
collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information Operations and Reports. 1215 Jefferson 
Davis Highway. Suite 1204, Arlington. VA 22202-4302, and to the Office of Management and Budget, Paperwork Reduction Project (0704-0188). Washington, DC 20503. 

1. AGENCY USE ONLY (Leave blank) 

2. REPORT DATE 

September 1994 

3. REPORT TYPE AND DATES COVERED 

Technical Memorandum 

4. TITLE AND SUBTITLE 

5. FUNDING NUMBERS 


A New Algorithm for Five-Hole Probe Calibration, Data Reduction, and 
Uncertainty Analysis 


6. AUTHOR(S) 

Bruce A. Reichert and Bruce J. Wendt 


WU- 505-62-52 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135-3191 


8. PERFORMING ORGANIZATION 
REPORT NUMBER 


E-8319 


9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Washington, D.C. 20546-0001 


10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 


NASA TM- 106458 


11. SUPPLEMENTARY NOTES 

Bruce A. Reichert, NASA Lewis Research Center; and Bruce J. Wendt, National Research Council-NASA Resident 
Research Associate at Lewis Research Center. Responsible person, Bruce A. Reichert, organization code 2660, 
(216) 433-8397. 


12a. DISTRIBUTION/AVAILABILITY STATEMENT 


Unclassified - Unlimited 
Subject Category 02 


12b. DISTRIBUTION CODE 


13. ABSTRACT (Maximum 200 words) 

A new algorithm for five-hole probe calibration and data reduction using a non-nulling method is developed. The 
significant features of the algorithm are: 1) two components of the unit vector in the flow direction replace pitch and yaw 
angles as flow direction variables, and 2) symmetry rules are developed that greatly simplify Taylor's series representa- 
tions of the calibration data. In data reduction, four pressure coefficients allow total pressure, static pressure, and flow 
direction to be calculated directly. The new algorithm's simplicity permits an analytical treatment of the propagation of 
uncertainty in five-hole probe measurement. The objectives of the uncertainty analysis are to quantify uncertainty of five- 
hole probe results (e.g., total pressure, static pressure, and flow direction) and determine the dependence of the result 
uncertainty on the uncertainty of all underlying experimental and calibration measurands. This study outlines a general 
procedure that other researchers may use to determine five-hole probe result uncertainty and provides guidance to 
improve measurement technique. The new algorithm is applied to calibrate and reduce data from a rake of five-hole 
probes. Here, ten individual probes are mounted on a single probe shaft and used simultaneously. Use of this probe is 
made practical by the simplicity afforded by this algorithm. 


14. SUBJECT TERMS 

How measurement; Pneumatic probes; Velocity measurement; Pressure measurement; 
Error analysis; Uncertainty analysis 


15. NUMBER OF PAGES 

20 


16. PRICE CODE 

A03 


17. SECURITY CLASSIFICATION 
OF REPORT 

Unclassified 


18. SECURITY CLASSIFICATION 
OF THIS PAGE 

Unclassified 


19. SECURITY CLASSIFICATION 
OF ABSTRACT 

Unclassified 


20. LIMITATION OF ABSTRACT 


NSN 7540-01-280-5500 


Standard Form 298 (Rev. 2-89) 
Prescribed by ANSI Std. Z39-18 
298-102 


National Aeronautics and 
Space Administration 

Lewis Research Center 

21000 Brookpark Rd. 
Cleveland, OH 44135-3191 


.! 




w « 

m > 

CD -r 
£ Q_ 


CQ — 
O <5 

5 = ® 

O Q. 


c 

i— 

0 

0 


CE 


o 

Q 


_o 

-Q 

2 

0 

> 

0 

T 3 

c 

Z> 


cc 

LU 

l— 

c o 
< 


CO 

o 

0L 


